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INTRODUCTION 


This  document  presents  a  summary  and  comparison  of  two  different  cross¬ 
correlation  techniques:  real  correlation  and  envelope  correlation  by  using  the  analytic 
signal.  Data  recorded  in  analog  form  are  subsequently  sampled  and  digitized  for  processing 
in  a  digital  computer.  The  correlator  operates  on  any  real  or  complex  signals  and  outputs 
the  magnitude,  phase,  real  part,  and  imaginary  parts  of  the  complex  output.  The  digitized 
data  can  be  directly  input  as  real  data  for  real  correlation,  or  they  can  be  input  in  the  form 
of  the  complex-valued  analytic  signal  for  an  envelope  correlation,  n 


REAL  CORRELATION 


The  data  are  input  into  the  correlator  with  only  real  parts,  as  shown  in  figure  1.  The 
correlator  assumes  the  imaginary  part  to  be  zero.  The  fast  Fourier  transform  (FFT)  of  each 
channel  is  performed  and  the  resulting  complex  spectra  are  Hermitian  (real  part  exhibit¬ 
ing  even  symmetry;  imaginary  part  exhibiting  odd  symmetry).  One  channel’s  spectrum  then 
undergoes  a  complex  conjugation.  Since  only  the  imaginary  parts  change  their  sign,  the 
result  is  still  Hermitian.  Then  the  two  spectra  undergo  a  complex  multiplication  and  the 
product  is  again  Hermitian.  Finally,  the  inverse  fast  Fourier  transform  (IFFT)  of  a  Hermi¬ 
tian  signal  returns  a  real-valued  sequence  of  data.  Since  the  output  is  completely  real,  there 
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Figure  1.  Real  correlation. 


ENVELOPE  CORRELATION 


Before  the  data  are  sent  to  the  correlator,  they  are  made  into  an  analytic  signal.  The 
analytic  signal  xa(t)  is  a  complex-valued  signal  formed  from  the  original  real  signal  x(t),  as 
shown  in  figure  2.  The  negative  frequency  values  of  the  spectrum  X(f)  are  set  to  zero  and 
the  positive  values  are  doubled.  We  have 

x(t)  -  X<f) 

xa(t)  -*2  -u(/)-X(J) 

where  u(J)  is  the  frequency  domain  unit  step  function.  Multiplication  in  the  frequency 
domain  results  in  convolution  in  the  time  domain  as 
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where  the  real  part  of  the  analytic  signal  is  just  x(r),  and  the  imaginary  part  x(t)  is  the 
Hilbert  transform  of  jc(r).  The  Hilbert  transform  pair  is  given  as 


x(t) 

t  -  T 


dr 


,  +00  V 

^-dr. 


The  analytic  signal  then  can  be  summarized  as: 


3t(x„)  =  x{t)  ,  3(.va)  =  x(t)  ,  |xa(/)|  =  yj  x2(t)  +  x2(0  , 


4>(r)  =  tan-1 


After  each  input  channel  has  been  converted  to  this  analytic  signal  representation, 
they  are  rent  through  the  correlator,  as  shown  in  figure  3.  The  complex-valued  spectra 
exhibit  no  symmetry  properties,  because  the  negative  region  of  frequency  is  zero-valued.  As 
a  result,  the  inverse  FFT  of  the  complex-conjugate  multiplication  is  complex-valued.  The 
resulting  correlation  has  phase  that  may  range  from  -ir  to  +tt.  The  magnitude  of  this  signal 


3 


ANALYTIC 


Figure  2.  Formation  of  the  analytic  signal 
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is  known  as  the  envelope  correlation.  It  is  also  interesting  to  note  that  the  real  part  of  it  is 
equivalent  to  the  ornput  of  the  real  correlation  method.  Therefore,  it  is  possible  to  obtain 
both  the  real  correlation  and  the  envelope  correlation  by  means  of  this  method. 


COMPARISONS 

It  is  instructive  to  further  examine  the  details  of  the  the  two  correlation  approaches. 
To  do  this,  we  consider  a  flat  rectangular  cross  spectrum  under  the  following  conditions: 

Case  1.  Correlation  of  nonanalytic  signals 

•  Full  band,  both  positive  and  negative  frequencies 

•  Bandpass,  both  positive  and  negative  frequencies. 

Case  2.  Correlation  of  analytic  signals. 

•  Full  band,  positive  frequencies  only 

•  Bandpass,  positive  frequencies  only. 

The  first  case  corresponds  to  real  correlation,  the  second  to  envelope  correlation. 
We  will  show  the  modulating  effects  that  bandpassed  signals  yield.  We  will  also  show  that 
we  can  obtain  the  real  correlation  by  the  methods  used  in  case  2. 

CASE  1:  REAL  CORRELATION  (NONANALYTIC  SIGNALS) 

In  general,  the  FFT  of  the  rectangular  pulse  is  a  sine  function,  as  shown  in  figures 
4a  and  4b. 


A  •  rect 


A  W  -  sine  ( Wt)  = 


AW-  sin  (uWt) 
nWt 


for  full  band. 


We  get  a  pure  sine  function  that  is  real-valued.  This  is  the  output  we  expect  where  the 
spectrum  is  flat  across  the  whole  band.  It  is  also  what  we  expect  for  basebanded  bandpass 
signal,  because  the  spectra  look  essentially  the  same.  We  now  examine  the  effects  of 
bandpassing  the  signal  (figures  4c  and  4d). 

6(f)  =  A  -  Rect  {^rj  +  A  '  Rect 

g(t)  =  AW  -  exp  ( j2irfct )  •  sine  ( Wt)  +  A  W  •  exp  (-j2nfct)  -  sine  ( Wt) 

=  AW  -  sine  ( Wt)[ex  p  ( j2irfct )  +  exp  (-y'27r/c/)] 

-2AW  -sin  (Wt)  -cos  (2wfet) 

2A  W  -sin  (irWt) 

= - — - •  cos  (2nfct)  for  bandpass. 

vWi 
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(c)  Cross  spectrum  (bandpass). 


(d)  Resulting  modulated  sine  function. 


Figure  4.  Real  correlation  (nonanalytic  signals). 


This  result  is  also  real-valued,  but  it  is  a  sine  function  that  is  modulated  by  a  cosine  of 
frequecy  fc.  We  can  see  this  modulation  and  the  effect  of  varying  W  and  fc  in  figures  5a 
through  e.  Whenever  W  -  2fc,  the  modulated  sine  becomes  a  pure  sine  with  twice  the 
frequency  of  the  envelope,  as  shown  in  figure  5f. 


CASE  2:  ENVELOPE  CORRELATION  (ANALYTIC  SIGNALS) 


In  general,  for  a  flat  rectangular  band  of  width  W  and  centered  at  fc  we  have 


G(f)  =  A  •  Rect 


f-fc 

w 


g(t)  -  A  tV  ■  exp  (j2irfct)  •  sine  ( Wt) 

-AW  sine  ( Wt)[  cos  ^2nfct)  +  j  •  sin  (27r/c/)] 

=  AW  -  [sine  ( We)  •  cos  (2  irfct)  +  j  ■  sine  ( Wt)  •  sin  (2irfct)]  • 

This  result  may  be  expressed  in  terms  of  magnitude,  phase,  real  part,  and  imaginary  part 
for  both  bandpass  and  full-band  (W  =  2 fc)  signals. 

|g(r)|  =  AW  [sine2  ( We)  ■  cos2  (27 rfct)  +  sine2  ( Wt)  •  sin2  (27r/c/)] 


=  AW-  sine  (Wt) 


for  full  band  and  bandpass. 
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=  tan"1  [tan  (271-/^)]  =  2irfct 


for  bandpass. 


W 

-2ir  —  -  it  Wt 
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for  full  band. 


=  AW  -  sine  ( Wt)  •  cos  (2irfct) 


for  bandpass. 


=  A  W  sin  (2TrWt) 
2t rWt 


-AW  sine  (2  We) 


for  full  band. 


-  AW-  sine  ( Wt)  •  sin  (2nfct) 

-  A  W  -sine  (Wt)-  sir.  (irWt) 


for  bandpass. 


for  full  band. 


These  full-band  and  bandpass  correlations  are  shown,  respectively,  in  figures  6  and  7.  In 
either  case,  the  magnitudes  are  equivalent  and  are  known  as  the  envelope  correlation.  It  is  a 
sine  function  dependent  only  on  the  bandwidth  W.  The  phase  of  the  correlation  shows 
rotations  through  27r  at  the  frequency  center  (fc)  of  the  band.  The  real  and  imaginary  parts 
are  sine  functions  modulated  by  sinusoidal  carriers  of  frequency  fc.  The  real  parts  are 
modulated  by  a  cosine,  and  the  Imaginary  parts  by  a  sine.  Since  we  are  dealing  here  with 
analytic  signals,  the  imaginary  part  is  just  the  Hilbert  transform  of  the  real  part.  As  we 
observed  earlier,  in  the  case  where  W  -  2fc  (full  band  or  basebanded  bandpass),  the  modu¬ 
lated  sine  becomes  a  pure  sine  of  twice  the  envelope  frequency.  It  is  also  interesting  to 
observe  that  the  real  parts  of  the  envelope  correlations  (case  2)  are  equivalent  to  the  corres¬ 
ponding  real  correlations  (case  1).  Therefore,  we  may  obtain  the  real  correlations’  results  by 
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means  of  the  envelope  correlation,  by  looking  at  the  real  part.  The  envelope  method  of 
processing  is  advantageous  in  that  it  gives  a  more  complete  and  thorough  representation  of 
the  correlation  function. 


EXAMPLES 

We  now  look  at  the  correlation  method,  using  actual  data  and  processing  software. 
First,  we  attempt  to  reproduce  the  results  obtained  previously  for  a  flat  power  spectrum. 
This  is  approximated  by  applying  software  (deline  and  delta  mean  equalizer)  to  flatten  the 
power  spectrum.  Then  the  data  are  bandpass  filtered  with  bandwidth  W  =  200  Hz  centered 
at  fc  -  400  Hz.  In  this  manner,  we  approximate  a  flat  rectangular  band  for  both  full  band 
(baseband)  and  bandpass,  as  shown  in  figure  8.  The  resulting  envelope  correlations’  mag¬ 
nitudes  and  phases  are  shown  in  figure  9.  We  can  also  see  the  real  and  imaginary  parts  in 
figure  10.  The  correlation  results  are  close  to  and  consistent  with  what  we  expect.  Figure  1 1 
shows  the  envelope  correlations  in  a  correlogram  display.  For  this  type  of  display,  the  min¬ 
imum  value  is  assigned  a  white  pixel,  and  the  maximum  value  a  black  pixel,  setting  the 
limits  of  the  gray  scale.  The  magnitude  correlograms  will  always  range  from  zero  to  1,  and 
the  phase  from  -v  to  +7r.  The  real  and  imaginary  correlograms  range  from  the  minimum 
correlation  (negative  values)  to  the  maximum  correlation  values. 

Next,  we  look  at  the  cross  correlation  between  signal  recorded  at  two  physically 
separated  locations.  We  view  them  in  terms  of  the  correlogram  outputs.  Figure  12  shows 
the  correlations  for  the  full  band,  performed  by  both  the  real  and  the  envelope  methods. 
We  observe  the  magnitude,  phase,  real,  and  imaginary  parts  of  the  correlation.  Again,  we 
see  that  the  real  part  of  the  envelope  correlation  is  equivalent  to  the  real  correlation.  The 
imaginary  part  of  the  real  correlation  is  zero,  with  any  minute  delations  spread  across  the 
entire  gray  scale.  Figure  13  shows  the  effects  of  bandpassing  the  signal.  We  see  both  the 
envelope  and  the  modulated  real  correlation.  Finally,  in  figure  14  we  see  the  effect  of  base¬ 
banding  the  input  data. 
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(e)  Basebanded. 
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Figure  8.  Approximation  of  flat  rectangular  spectrum,  using  actual  data. 


(c)  Phase  of  bandpass. 


(d)  Phase  of  baseband. 

Figure  9.  Autocorrelation  of  figure  8  (1). 
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(c)  Imaginary  part  of  bandpass. 


Figure  10.  Autocorrelation  of  figure  8  (2). 
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Figure  11.  Envelope  correlograms  of  figure  8  (bandpas,  and  basebanded). 
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Figure  14.  Sensor  data  correlograms,  basebanded  (real  and  envelope) 


SUMMARY 


We  have  described  two  methods  of  processing,  real  and  envelope,  that  can  be  used 
to  correlate  signals.  It  is  possible  to  obtain  the  real  correlation  with  the  envelope  correla¬ 
tion  method.  It  is  advantageous  to  process  with  the  analytic  signal  because  it  p  .  vides  both 
of  these  correlation  functions.  With  both  correlation  functions  available,  one  may  then 
choose  which  is  better  suited  for  the  job. 
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